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PREVIOUS  WORK 


The  determination  of  the  coordinates  (X  and  Y)  of  the  earth's  spin  axis 
from  Dqppler  observation  has  been  described  by  Anderle  and  his  colleagues 
in  a series  of  publications  and  reports. 

The  method  of  computation  was  briefly  explained  by  Anderle  and  Beuglass 
(1970).  A more  detailed  description  of  the  observational  procedures,  the  data 
reduction  techniques,  and  error  sources  was  given  by  Anderle  (1973a) . 

Results  of  Doppler  data  analysis  and  comparison  with  other  determinations 
are  discussed  in  numerous  places.  All  Doppler  results  are  based  on  five-day 
mean  values  of  X and  Y.  They  are  tabulated  and  discussed  for  1969  by  Anderle 
and  Beuglass  (1970);  for  1967  to  1970  by  Anderle  (1970);  for  1969  to  1971  by 
Anderle  (1972) ; for  1972  by  Anderle  (1973b) ; and  for  1973  by  Beuglass  (1974) . 
The  five-day  means  for  the  years  1974  to  1977  are  given  in  this  report.  A few 
tv«-day  solutions  for  1964  to  1969  may  be  found  in  Anderle  (1973b) , Appendix  F. 

The  above  data  are  normally  shown  as  plots  of  X versus  time,  Y versus 
time,  aro  X versus  Y.  They  are  given  by  Anderle  and  Beuglass  (1970)  for  1969; 
by  Anderle  (1970)  for  1967,  1968,  1969  and  1970;  by  Beuglass  and  Anderle  (1972) 
for  1970;  by  Anderle  (1972)  for  1969,  1970,  and  1971;  by  Anderle  (1973b)  for 
the  period  1964  through  1967  and  for  1972;  by  Anderle  (1973a)  for  mid-1971 
to  mid-1972;  by  Beuglass  (1974)  for  1973;  and  by  Anderle  {1976b)  for  1975. 

Plots  for  the  years  1974  to  1977  are  shown  in  this  report. 

Anderle  (1976a)  has  also  compared  Doppler  derived  pole  coordinates  with 
classical  optical  solutions.  He  plots  the  differences  BIH-ILS,*  DMA  (Doppler) 
-BIH,**  and  DMA-ILS  for  the  span  1964  to  1975.  He  also  tabulates  yearly  mean 
values  for  above  differences  as  well  as  associated  statistics.  Anderle  (1976b) 
adds  the  comparison  DMA-IPMSt  and  shows  more  detail  by  breaking  the  plots  into 
two  spans,  1964  through  1969  and  1970  to  1975. 


OBSERVATIONS  AND  DATA  REDUCTIONS 


The  following  is  a very  brief  description  of  observational  data  and  their 
analyses.  Details  may  be  found  in  the  references  listed  at  the  end  of  this 
report,  especially  in  Anderle  (i973a)  and  (1976b) . 

Observations  are  the  Doppler  shifts  in  the  continuous  radio  frequencies 
at  150  and  400  MHz  transmitted  by  the  U.  S.  Navigation  System  satellites 


* BIH-ILS  is  Bureau  International  de  L'Heure-lnternational  Latitude  Service 

**  DMA-BIH  is  Defense  Mapping  Agency-Bureau  International  de  L'Heure 
t DMA-IPMS  is  Defense  Mapping  Agency-International  Polar  Motion  Service 
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(Kershnec,  1967).  Analog  combination  of  these  two  frequencies  permits  elimi- 
nation of  one  large  error  source;  namely,  the  first-order  ionospheric  refrac- 
tion effect. 

The  number  of  satellites  being  observed  varies  between  two  and  five, 
depending  on  Navy  requirements.  Table  1 shows  which  satellites  were  observed, 
and  when,  for  the  years  1974  to  1977. 


Table  1.  Available  Dc^pler  Satellite  Data 
(Day  Numbers) 


Satellites 


1967-34A 

1967-48A 

1967-92A 

1970-67A 

1973-81A 

1974 

166-280 

1-87 

89-363 

1975 

2-362 

13-363 

1976 

155-365 

6-364 

1-157 

1977 

7-167 

21-167 

21-167 

6-364 

21-365 

Observations  are  taken  by  as  many  ground  stations  as  are  operational. 

They  increased  in  number  from  about  13  in  1969  (Anderle  and  Beuglass,  1970) 
to  about  20  in  recent  years  (Anderle,  1976b) . 

All  observations  taken  during  a 48-hr  time  span  are  used  in  a least-squares 
solution  to  improve,  primarily,  the  orbital  parameters.  During  this  process, 
the  satellite  orbits  are  numerically  integrated  by  Cowell's  method,  that  is, 
the  Gauss-Jackson  algorithm  applied  to  the  differential  equations  in  the  rect- 
angular accelerations.  The  program  is  normally  run  with  a 60-sec  integration 
step  size  and  order  12.  The  reference  frame  is  the  mean  equator  and  equinox 
at  the  beginning  of  the  observation  span.  The  mathematical  model  contains 
about  480  gravity  terms,  atmospheric  drag,  radiation  pressure,  luni-solar  solid 
earth  tides,  with  the  Love  number  presently  set  at  0.26.  The  force  field  is 
complete  enough  to  permit  determination  of  the  satellite's  position  good  to 
about  1 m. 

The  solution  also  contains,  among  other  parameters,  the  coordinates  X and  Y 
of  the  spin  axis,  referred  to  the  Conventional  International  Origin  (CIO) . 

Such  two-day  solutions  are  obtained  separately  for  each  satellite.  Subsequent- 
ly, all  two-day  solutions  from  up  to  five  different  satellites  are  combined  into 
five-day  means.  The  latter  is  published  by  the  U.  S.  Naval  Observatory  in 
"Preliminary  Times  and  Coordinates  of  the  Pole,  Series  7." 

The  computation  of  pole  positions  based  on  Doppler  observations  originated 
at  the  Naval  Weapons  Laboratory  (now  the  Naval  Surface  Weapons  Center  (NSWC) ) . 

In  April  1975,  the  responsibility  of  computing  Navy  navigational  satellite 
(NAVSAT)  orbit,  and,  hence,  the  derivation  of  pole  positions  was  transferred 
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to  the  Topographic  Center  of  the  Defense  Mapping  Agency  (EMATC) . Since  EMATC 
employs  the  same  computer  programs,  the  transfer  did  not  affect  position  results. 

Over  the  years,  there  have  been  a number  of  changes  in  the  observation 
station  network  and  observation  techniques  (Anderle,  1973a)  as  well  as  improve- 
ments in  the  data  reduction  methods  (Anderle,  1972).  However,  the  procedures 
have  been  essentially  the  same  since  August,  1971,  so  that  Doppler  results 
after  this  date  are  believed  to  be  homogeneous. 


DOPPLER  POLE  POSITIONING  ACCURACY 


The  formal  standard  deviation  for  the  polar  coordinates  from  a two-day 
solution  is  about  5 cm  during  the  second  half  of  1977.  But  it  must  be  re- 
membered that  such  solutions  are  made  for  each  satellite  separately.  All  two- 
day  solutions  are  then  combined  into  five-day  means.  Subsequently,  one  can 
compute  the  more  realistic  standard  deviation  of  a two-day  coordinate  with 
respect  to  the  five-day  mean.  That  number  is  presently  a bit  less  than  40  cm. 
The  standard  deviation  of  the  five-day  mean  itself  (standard  error)  has  been 
just  under  20  cm  for  the  last  two  years. 

The  increase  in  accuracy  from  1967  to  1977  is  shown  in  Table  2.  However, 
the  data  before  and  after  1972  are  not  immediately  comparable.  Polar  co- 
ordinates until  August  1971  were  obtained  using  the  tangential  c(»nponent  of 
the  station  navigations  only,  and  they  were  one-day  solutions.  Moreover,  they 
were  computed  after  orbit  improvement,  not  in  a simultaneous  least-squares 
solution. 

Table  2.  Preliminary  Yearly  FMS  of  Standard  Deviations  and  Errors 
Standard  Deviations  (m)  Standard  Errors  (m) 


(two-day  solutions*) 

(five-day  means) 

X 

Y 

Av. 

X 

Y 

Av. 

1967 

1.65 

1.78 

1.72 

0.89 

0.74 

0.82 

1968 

1.48 

1.60 

1.54 

0.86 

0.93 

0.90 

1969 

1.51 

1.27 

1.40 

0.69 

0.60 

0.65 

1970 

1.25 

1.15 

1.20 

0.57 

0.53 

0.55 

1971 

1.16 

1.39 

1.28 

0.52 

0.62 

0.57 

1972 

0.75 

0.69 

0.72 

0.37 

0.32 

0.35 

1973 

0.38 

0.44 

0.41 

0.22 

0.28 

0.25 

1974 

0.48 

0.50 

0.49 

0.30 

0.32 

0.31 

1975 

0.47 

0.36 

0.42 

0.22 

0.18 

0.20 

1976 

0.40 

0.30 

0.35 

0.20 

0.15 

0.18 

1977 

0.41 

0.34 

0.38 

0.18 

0.15 

0.17 

*One-day  solutions  before  1972 
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Anderle  (1973a)  pointed  out  that  the  principal  error  source  in  Doppler 
polar  coordinates  is  due  to  inadequate  knowledge  of  the  earth  gravity  field. 
Despite  recent  advances,  this  remains  true  today. 


RESULTS  1974-1977 

Anderle  and  his  colleagues  have  already  published  diagrams  and  tables 
summarizing  polar  motion  during  1974  and  1975.  Since  some  of  their  results 
were  based  on  preliminary  data,  they  ate  repeated  here  using  final  data. 
Final  values  were  also  available  for  1976  while  some  1977  results  are  still 
preliminary.  They  will  be  identified  as  such  below. 


TABULATION  OF  DOPPLER  RESULTS 


Table  3 is  a sample  containing  the  two-  and  five-day  Doppler  solutions 
for  polar  coordinates.  The  ccmplete  tables  for  the  years  1974  to  1977  will  be 
published  in  a forthcoming  NSWC  report.  The  first  two  columns  show  the  day 
numbers  for  each  two-day  solution.  They  are  followed  by  X and  Y and  their  formal 
standard  deviations  (labelled  "Standard  Error")  as  obtained  from  the  covari- 
ance matrix  of  the  least-squares  solution.  The  last  two  columns  are  the 
satellite  designation  and  the  nominal  value  for  UTC-UTl.  The  latter  informa- 
tion is  not  used  in  our  pole-position  calculations. 

In  the  last  three  lines  of  each  block,  only  the  first  three  columns  are  of 
interest.  The  first  column  shows  the  day  number  for  which  the  five-day  average 
is  being  competed.  Columns  two  and  three  show  the  weighted  averages  for  X and 
Y,  where  the  weight  is  taken  as  l/o^,  o being  the  two-day  standard  deviations 
mentioned  above.  The  line  marked  STD  DEV  is  the  weighted  standard  deviation 
of  a two-day  solution  with  respect  to  the  five-day  mean.  The  last  line,  labelled 
STD  ERR,  is  the  previous  line  divided  by  the  square  root  of  n.  It  is,  therefore, 
the  standard  deviation  of  the  five-day  mean. 

Note  that  the  program  is  presently  limited  to  include  only  the  first  four 
two-day  solutions  for  any  given  day  in  the  five-day  means,  even  though  ail 
available  two-day  results  are  listed. 


Table  3.  Dahlgren  Polar  Monitoring  Service 
NWL  9 Pole  Report  Revision 


Dally  Solution  (m) 


Bi-Daily  (m) 


Pole  Position 

Standard 

Error 

Solution 

Nominal 

Days 

1977 

X 

Y 

X 

Y 

X Y 

Satellite 

OTO-UTl 

333. 

334. 

2.58 

1.11 

.05 

.05 

.01-947058.66 

1973  81A 

244000.00 

334. 

335. 

3.46 

2.48 

.04 

.05 

.01-953246.60 

1970-67A 

247000.00 

335. 

336. 

4.21 

1.12 

.05 

.05 

.01-947058.66 

1973  SLA 

250000.00 

336. 

337. 

3.00 

2.30 

.05 

.06 

.01-953246.60 

1970-67A 

254000.00 

337. 

338. 

2.72 

1.74 

.06 

.05 

.01-947058.66 

1973  8LA 

258000.00 

HEAN 

336. 

3.25 

1.75 

.05 

.05 

STD  DEV 

336. 

.66 

.65 

.00 

.00 

STD  ERR 

336. 

.30 

.29 

338. 

339. 

2.74 

1.18 

.05 

.05 

.01-953246.60 

1970-67A 

262000.00 

339. 

340. 

2.58 

1.44 

.05 

.05 

.01-947058.66 

1973  61A 

266000.00 

340. 

341. 

2.61 

1.24 

.05 

.05 

.01-953246.60 

1970-67A 

270000.00 

341. 

342. 

2.17 

1.91 

.05 

.05 

.01-947058.66 

1973  81A 

274000.00 

342. 

343. 

3.34 

1.06 

.05 

.05 

.01-953246.60 

1970-67A 

276000.00 

MEAN 

341. 

2.69 

1.  J6 

.05 

.05 

STD  DE' 

Ml. 

.41 

.33 

.00 

.00 

STD  ERl 

341. 

.18 

.15 

34  3. 

344. 

1.63 

.45 

.06 

.05 

.01-947058.66 

1973  81A 

279000.00 

344. 

345. 

1.94 

.52 

.07 

.06 

.01-953246.60 

1970-67A 

282000.00 

345. 

346. 

2.06 

.46 

.05 

.06 

.01-947058.66 

1973  81A 

285000.00 

346. 

347. 

2.05 

1.38 

.05 

.05 

.01-953246.60 

1970-67A 

288000.00 

347. 

348. 

1.26 

.51 

.06 

.05 

.01-947058.66 

1973  81A 

290000.00 

MEAN 

346. 

1.82 

.66 

.06 

.05 

STD  DEV 

346. 

.34 

.40 

.01 

.00 

STD  ERR 

346. 

.15 

.18 

348. 

349. 

2.29 

1.99 

.05 

.05 

.01-953246.60 

1970-67A 

294000.00 

349. 

350. 

1.40 

.58 

.05 

.05 

.01-947058.66 

1973  81A 

295000.00 

350. 

351. 

.81 

.83 

.05 

.06 

.01-953246.60 

1970-67A 

297000.00 

351. 

352. 

1.17 

.77 

.05 

.05 

.01-947058.66 

1973  81A 

299000.00 

352. 

353. 

1.52 

.49 

.05 

.05 

.01-953246.60 

1970-67A 

300000.00 

MEAN 

351. 

1.45 

.89 

.05 

.05 

STD  DEV 

351. 

.56 

.59 

.00 

.00 

STD  ERR 

351. 

.25 

.27 

POLE  COORDINATE  PLOTS 


The  motion  of  the  pole  during  the  years  1974  to  1977  may  be  seen  at  a glance 
in  Figures  1 through  4.  The  most  striking  feature  is  the  increase  in  amplitude, 
from  about  3 to  8 m (half  amplitude) . It  is  simply  a consequence  of  superposition 
of  the  365.25-  and  432-day  components,  out  of  phase  in  1974,  and  in  phase  in  1977. 

The  Doppler  data,  now  labelled  DMA,  are  easily  identified  by  their  1 
error  ellipses.  These  are  the  STD  ERR  of  the  five-day  means  shown  in  Table  3. 
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Also  shown  arc  the  polar  coordinates  frcxn  three  other  sources,  namely 
BIH,  ILS,  and  IPMS.  They  are  plotted  as  solid  lines,  dashes,  and  alternating 
dots  and  dashes,  respectively.  As  in  earlier  years,  the  agreement  between 
BIH  and  DMA  is  quite  good.  However,  it  must  be  pointed  out  that  the  Doppler 
data  are  used,  in  addition  to  optical  observations,  in  deriving  the  BIH  pole 
position  results  quoted  for  1977.  The  agreement  between  IPMS  and  BIH  or 
Doppler  is  reasonably  good;  the  difference  is  2 m only  early  in  1976,  but 
usually  much  less.  ILS,  however,  frequently  differs  from  the  other  three 
determinations.  The  discrepancy  reaches  3.5  m at  several  times.  The  total 
excursion  of  the  pole,  as  shown  by  ILS,  appears  to  be  somewhat  smaller  than 
for  the  other  three  determinations. 

The  BIH  path  shown  in  Figure  4 appears  exceedingly  smooth.  This  is  due 
to  the  fact  that  only  "smoothed"  data  were  available  at  the  time  the  above 
plots  were  prepared.  In  the  meantime,  the  1977  Circular  D "raw"  data  were  ac- 
quired and  used  in  Figure  5.  It  is  believed  to  be  more  representative  of  the 
accuracy  of  the  BIH  five-day  values. 

Figures  6 to  13  permit  a comparison  of  the  various  polar  motion  services 
separated  into  the  X and  Y coordinates.  It  may  be  seen  that  the  differences 
are  always  larger  in  X than  they  are  in  Y.  In  fact,  for  1974  through  1976  the 
agreement  of  all  four  services  in  Y is  remarkably  good.  Note  that  again  the 
smoothed  BIH  values  were  used  in  plotting  Figures  9 and  13. 


Figure  1.  1974  Pole  Path 
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Figure  2.  1975  Pole  Path 
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Figure  12.  Pole  Position  for  1976,  Y Ccroponent 


Figure  13.  Pole  Position  for  1977,  Y Component 
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DIFFERENCES  IN  POLE  COORDINATES 


Anderle  (1976b)  published  plots  of  differences  in  the  X and  Y coordinates 
of  the  pole  for  the  time  spans  1964  through  1969  and  1970  through  1975. 

Similar  plots,  Figures  14  and  15,  are  given  in  this  report  for  the  interval 
1974  through  1977.  These  diagrams  show  the  differences  in  the  four  pairs 
BIH-ILS,  DMA-ILS,  DMA-BIH,  and  DMA-IPMS  quite  clearly.  By  and  large,  the 
Y-coordinates  agree  well,  except  for  the  ILS  excursions  in  1976  and  1977. 

In  X,  however,  all  four  pairs  show  significant  biases.  ILS  again  shows  some 
large  variations  with  respect  to  BIH  and  DMA. 

Tables  4 and  5 are  a continuation  of  similar  information  published  by 
Anderle  in  earlier  reports.  They  list  the  yearly  average  difference  for  each 
of  the  four  pairs  being  compared,  as  well  as  the  standard  deviation  of  the 
individual  difference  with  respect  to  the  annual  mean.  Individual  points  in- 
volving either  ILS  or  IPMS  would  be  18  days  apart,  while  DMA-BIH  is  formed 
every  five  days.  Footnotes  to  Tables  4 and  5 contain  additional  information 
concerning  data  sources  and  reference  frames. 


THE  CHANDLER  PERIOD 


It  is  well  known  that  the  principal  periodic  contents  of  the  motion  of 
the  pole  are  the  Chandler  period  and  the  annual  term.  In  order  to  determine 
the  former,  Anderle  (1977)*  adapted  an  existing  program  to  fit  to  the  data  an 
expression  of  the  form 


and  a similar  equation  for  Y.  P^  is  the  unknown  Chandler  period,  and  the  X^, 
A,  and  C are  five  numerical  coefficients  to  be  determined  by  least-squares 
fits.  One  assumes  a value  for  P^,  obtains  am  expression  for  X,.f>n,p  and  Y^g^p, 
and  forms  the  residuals  and  their  BMS.  This  is  repeated  for  sevec^al  values^ 
of  Pq,  and  a parabola  is  fitted  to  three  such  pairs  of  points.  Finally,  one 
computes  the  value  of  P^  for  which  the  rms  parabola  has  its  minimum.  Obvious- 
ly, to  obtain  P^  directly  from  a least-squares  solution  is  more  elegant,  but 
the  above  procedure  permitted  the  use  of  existing  coding. 


* Anderle,  R.  J.,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory,  Dahlgren, 
Virginia,  office  memo,  1977. 
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Mean  Differenc«  C) Standard  Deviation  of  Difference  C) 

X-Coordinate  Y-Coordinate  X-Coordinate  Y-Coordinate 


(5)  Mean  corrections  of  -2.06  and  1.48  m were  added  to  1967  NW1.8D  X-  and  Y-pole  positions, 
respectively,  based  on  conparisons  with  12  NHL90  results  given  in  TR-2734  and  2952. 


Table  5.  Coordinate  Systems  and  Gravity  Fields 


YEAR 

COORDINATE  SYSTEM 

GRAVITY  FIELD 

DMA  POLE  POSITION 

1964 

NWL-9D 

NWL-9B 

NHL  TR-2734, 

2952  (9-Mean  Positions) 

1965 

NWL-9D 

NWL-9B 

NHL  TR-2734, 

2952  (13-Hean  Positions) 

1966 

NWL-9D 

NWL-9B 

NHL  TR-2734, 

2952  (6-Mean  Positions) 

1967 

NWL-8D*^> 

NVn,-8B 

NWL-8D  (20  Feb)  (5) 

Preprint 

1968 

NWL-8F  (19  Jan) 

NVn,-9H  (18  Apr) 

Preprint 

1969 

NWI,-8F*^'®I 

NWL-'8H<'^'®> 

NHL  TR-2734 

1970 

NWL-9C  (20  Dec) 

NWL-9B  (13  Feb) 

NHL  TR-2734 

1971 

NWL-9D  (18  Oct) 

NWL-9b'^' 

NHL  TR-2734 

1972 

NWL-9D 

NWL-9B 

NHL  TR-2952 

1973 

NWL-9D 

NVfL-lOE  (2  Jan) 

NHL  TR-3181 

1974 

NWL-9D 

NWL-lOE 

DMA  Heelcly  Reports 

1975 

NWL-9D 

NWL-lOE 

DMA  Weekly  Reports 

1976 

NWI.-9D 

NHL- 10  E 

DMA  Weekly  Reports 

1977 

NWI,-9Z-2  (15  June) 

NWL-lOE 

DMA  Weekly  Reports 

(6)  Mean  corrections  of  -2.37  and  2.04  m were  added  to  1968  NWL9D  X-  and  Y-pole  positions, 

respectively,  based  on  comparisons  with  11  NHL9D  results  given  in  TR-2734  and  2952. 

(7)  Mean  corrections  of  -0.07  and  -2.35  m were  added  to  1969  NWL8D  X-  and  Y-pole  positions, 

respectively,  based  on  comparisons  with  12  NHL9D  results  given  in  TR-2734  and  2952. 

(8)  TR-2734  gives  pole  positions  for  1969-1970  computed  after  adjusting  NVn,8F  latitude  residuals 

to  NNLIOD  system 

(9)  DHA  pole  positions  for  1964-1966  and  after  August  1971  are  based  on  simultaneous  solution 
for  orbit  constants  and  pole  position  rather  than  sequential  solutions. 
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Figures  16  and  17  depict  the  curves  discussed  above.  The  legend  shows 
that  12  years  of  data  were  used  for  the  three  astronomical  sources/  while  only 
six  years  of  Doppler  data  were  available  for  this  analysis.  It  can  be  seen 
that  the  astronomical  services  agree  quite  well.  The  minima  are  near  432  days, 
and  they  are  well-defined.  The  Doppler  X-coordinate  curve  is  quite  useless, 
and  that  for  the  Y-coordinate  is  of  dubious  value.  It  was  quickly  found  that 
the  short  six-year  time  span  is  responsible.  Solutions  for  the  three  astronomi 
cal  sources  over  six  years  produced  results  comparable  to  the  Doppler  curves. 


«r0  4/5  4S0  433  440 

ASSunKO  CHMutOLitt  ftmoo  m O^rs 


Figure  16.  Residuals  after  S-Parameter  Fit  X-Ccmponent 
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Figure  17.  Residuals  after  5-Paranieter  Fit  Y-Conponent 


Table  6 contains  the  results  of  the  P^,  computations  explained  above.  In 
obtaining  the  averages  and  the  mean  valuer  unit  weight  was  assumed.  The  error 
bound  of  0.2  day  was  calculated  from  the  scatter  of  the  six  individual  values. 


Table  6.  Chandler  Period  * 
Pj,  Days 


(X) 


ILS  431.71 
BIH  431.94 
IPMS  431.77 


(Y) 


432.87  432.29 

432.00  431.97 

431.79  431.78 


* Mean  Value:  P^  = 432.0  0.2  days 

Although  of  questionable  value,  P^  was  also  computed  from  the  DMA  Y-curve. 
It  yields  432.2  days,  in  reasonable  agreement  with  our  adopted  values  of  432.0 
days. 


Our  determination  is  also  in  good  agreement  with  Markowitz  (1976) , who 
obtains  432.02  0.15  days.  It  compares  reasoneOsly  well  with  Vicente  and  Currie 

(1976),  who  quote  433.2  0.8  days. 

The  above  mentioned  residuals  are  believed  to  contain  other  periodicities. 
Bowman  and  Leroy  (1976)  , among  others,  have  performed  a spectral  analysis  of 
the  X-  and  Y-components  themselves,  with  the  following  results: 


Table  7.  Bowman/Leroy  Spectral  Analysis 


Frequency 

(cycles/year) 


Period  Amplitude 
(days)  (m) 


0.85 

430 

5.84 

0.6 

1.0 

365 

4.84 

0.6 

1.3 

280 

0.49 

0.6 

2.0 

180 

0.23 

0.6 

2.5 

145 

0.12 

0.6 

4.0 

90 

0.11 

0.6 

Their  analysis  is  b2sed  on  five  years  of  Doppler  data.  Considering  our 
earlier  difficulties  with  such  a short  time  span,  perhaps  considerable  strength 
could  be  added  to  the  soljtion  by  including  astronomical  data.  A particularly 
attractive  time  span  would  be  13  years,  corresponding  to  almost  exactly  11 
Chandler  cycles.  However,  reliable  Dqppler  data  does  no:  yet  exist  for  such 
an  interval. 


ADVANTAGES  AND  DI SADV7UJTAGES  OF  DOPPLER 


Doppler  observations  are  taken  day  and  night,  and  under  any  cloud  covet. 
This  all-weather  capability  is  one  of  its  major  assets.  Doppler  data  are 
also  less  sensitive  to  tropospheric  effects  than  are  optical  observations. 
Moreover,  they  are  independent  of  star  catalog  position  errors.  Perhaps 
Doppler's  greatest  value  lies  in  the  fact  that  it  adds  a totally  independent 
pole  position  determination  to  the  classical  methods. 

Systematic  errors  due  to  an  inadequate  knowledge  of  the  gravity  field  is 
the  major  disadvantage  of  Doppler.  Results  are  also  affected  by  changes  in 
the  station  network  and  atmospheric  drag  variations  during  a two-day  span. 
Computing  Doppler  pole  positions  is  quite  expensive.  At  the  present  time, 
however,  they  are  obtained  as  by-products  in  orbit  improvement  runs  performed 
by  DMA.  Finally,  although  TRANSIT  satellites  have  shown  remarkable  endurance, 
their  life  time  is  finite. 


FUTURE  PLANS 

The  planning  of  drag-free  satellites  is  underway  at  the  Applied  Physics 
Laboratory,  Johns  Hopkins  University.  Once  operational,  effects  due  to  drag 
would  be  eliminated  and  a better  gravity  field  could  be  determined,  resulting 
indirectly  in  better  orbits  and  pole  positions. 

The  earth  gravity  field  is  continuously  being  improved,  especially  by  the 
National  Aeronautics  and  Space  Administration  (NASA) . NSWC  has  also  begun  work 
on  a major  new  geodetic  solution.  Other  improvements  in  the  mathematical  model 
are  planned,  especially  better  representations  of  the  various  tide  effects. 


SUmARY 


Computations  of  polar  coordinates  from  Doppler  observations  have  been  per- 
formed in  recent  years  by  DMA.  During  the  first  half  of  1977  as  many  as  five 
satellites  were  observed.  The  standard  deviation  of  a cwo-day  polar  coordinate 
solution  is  now  better  than  40  cm,  that  for  the  five-day  mean  under  20  cm. 
Agreement  between  the  four  services  ranges  from  excellent  to  only  fair.  There 
are  no  significant  problems  in  the  Y-coordinate,  except  a 1.5  m standard 
deviation  in  1977  for  comparisons  involving  ILS.  The  X-coordinate  shows  both 
large  biases  emd  standard  deviations. 

It  is  found  that  six  years  of  Doppler  data  are  not  enough  to  derive  a 
reliable  Ch^uldler  period.  Hence,  12  years  of  data  from  the  three  astronomical 
services  were  taken  to  compute  a Chandler  period  of  432.0  0.2  days.  Residuals 

suggest  the  existence  of  additional  periodic  terms. 
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